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Abstract 

■ We investigate the Hamiltonian analysis of Nappi-Witten model (WZW action 

based on non semi simple gauge group) and find a time dependent non-commutativity 
by canonical quantization. Our procedure is based on constraint analysis of the model 
in two parts. A first class analysis is used for gauge fixing the original model following 
by a second class analysis in which the boundary condition are treated as Dirac 
constraints. We find the reduced phase space by imposing our second class constraints 
<]J ' on the variables in an extended Fourier space. 
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£j 1 Introduction 

Treating boundary conditions as Dirac constrains has been considered in the recent years 
by so many authors [TJ |2j [3j H] . This approach has been used first in studying the Polyakov 
string coupled to a B-field. The common feature of all works is non commutativity of the 
coordinate fields on the boundaries which may lie on some brains, as first predicted by 
[Sj. However, there are different approaches in defining the constraints and investigating 
their consistency in time. We have reviewed the whole subject in our previous work [6] 
and showed if we impose the set of constraints on the Fourier expansions of the fields, the 
redundant modes will be omitted in a natural way. 

For simple physical models obeying linear equations of motion, the ordinary Fourier 
expansion gives appropriate coordinates to reach the reduced phase space. In other words, 
the infinite set of second class constraints emerging as the result of boundary conditions, 
forces us to omit a number of Fourier modes. However, ordinary Fourier transformation 
is not essential for quantization; it is just one tool that acts well for most physical models 
at hand. In the general case one should search for "appropriate coordinates", in which 
imposing the set of second class constraints is equivalent to omitting some canonical pairs 
from the theory. 

In this paper we study the constraint structure of the Nappi-Witten model in the Hamil- 
tonian formalism. This model acquires complicated boundary conditions so that the or- 
dinary Fourier expansion seems inadequate to impose the whole set of constraints which 
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emerge from the boundary conditions. Nevertheless, the Nappi-Witten model, on its own 
grands, is an attractive one since it describes a non semi simple gauge group as well as giving 
time dependent non commutativity in some gauges [7j. Our next interest is to emphasize 
that solving the equations of motion is not necessarily needed for quantizing a theory; the 
only necessity is finding the dynamics of the constraints and construct their algebra with 
the Hamiltonian such that they remain consistent with time on the constraint surface. 

We give a precise Hamiltonian treatment of the model in which the constraint structure 
is followed step by step from the initial action to the final reduced phase space. In section 2 
we introduce the model and find primary and secondary constraints of the system. Section 
3 is devoted to fixing the gauge by introducing appropriate gauge fixing conditions. In 
section 4 we follow our strategy of treating the boundary conditions as primary Dirac 
constraints and follow their consistencies. The boundary conditions which come from the 
original action, in fact, make the system more complicated. So, it is not possible to write 
down the solutions in a closed form similar to a simple Fourier expansion (see reference 
[8]). We try to find a basis which is appropriate for imposing the infinite set of constraints 
in section 5. In section 6 we will give our concluding remarks and will compare our results 
with parallel approaches. 



2 Hamiltonian structure of the model 

The Nappi-Witten model describes a 4-component bosonic string X a = (ai,a 2 ,u,v) living 
in the background metric G a b(X) and coupled to a S-field. The action is given as: 



5 



d 2 a 



-gg^G^d^djX' + B ab e ij d^djX 1 



(1) 



where 



G(X) 



( 1 





02 
2 


o\ 





1 


ai 
2 





12 
2 


ai 
2 


b 


1 







1 


v ) 



,B(X) 



( 


a 





o \ 


—u 































o J 



(2) 



The special form of G(X) and B(X) are chosen so that the gauge group of the model is 
non semi-simple [8]. The metric field can be written in terms of the following variables: 



^2 



r, ^3 



v / (ff 01)2_ p 00 ff ll' 



In terms of the variables X a and N Q the action becomes: 
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— G ab (X)(X a X b - 2N 2 X a X' b + (Nl - Nl)X' a X' b ) + 2B ab X a X' b 
Ni 



(3) 



(4) 



where dot and prime means temporal and spatial derivatives, respectively. The canonical 
momenta 7i a and p a conjugate to N a and X a are: 



n a = a = 1,2,3 

Pi = ^(2a; + ueijdj) - ^(2a- + u'e^aj) + 2ue ij a' j 

l.*! Mr, ^ / 

Pv 



p u = ^(2bii + 2v + eijCtidj) - j^{2bu' + 2v' + e ij a' i a j/ 



(5) 
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The Canonical Hamiltonian reads: 



H= j d 2 a±-G ab (F a F b - (iV 2 2 - N 2 )X ,a X ,b ), (6) 

where 

F a = X a = Ni(G~ 1 ) ab (p b - B bc X' c ) + N 2 B ab X' b (7) 
In terms of component fields dj, u and v we have 

H = J d 2 a{N^> 1 + iV 2 tt 2 ) (8) 

where 

^ = \v\ + \tijVvaiPj + \VuVv -\bp 2 v + jQa 2 pl 

+u / €ija' i aj + SijUa'tPj + \up v d { ai + (1 + u 2 )af + bu' 2 + 2u'v' (9) 
tt 2 = a'^ + u'p u + v'p v , 

As can be seen from Eqs. (JSJ) the momenta n a are primary constraints. The dynamics of 
the system is achieved by the total Hamiltonian: 

H T = H + J da\ a n a (a, r), (10) 

in which X a are Lagrange multipliers. As usual we should impose the consistency conditions 
on the constraints so that they remain valid during the time. For this reason we demand 
7r Q m 0, where ~ means weak equality i.e. equality on the constraint surface. Using Eqs. 
( [10]) and fl6]) we have: 

tt 1 = {ir\H T } = -tt 1 

p = {n 2 ,H T } = -tt 2 (11) 
tt 3 = {tt 3 , # T } = 0, 

Therefore, the consistency of three primary constraints ir a gives two second level constraints 
tt 1 and tt 2 . In this way we have so far two levels of constraints as 

In order to investigate the consistency of second level constraints, we need to calculate the 
Poisson brackets of VI/ 1 (cr, r) and tt 2 (a, r) at different points. Direct calculation, using the 
fundamental Poisson brackets among the four conjugate pairs (u,p u ), (v,p v ) and (a;,pj) 
gives: 

{tt >, t), tt V, r)} = ±(tt 2 (a, r)d a - tt 2 (a', r)d al )5{a - a') 

V, r ), V 2 (a', t)} = V, r)dj(a - a') (13) 
{tt 2 (a, r), tt 2 (a', r)} = ±(tt 2 (a, r)d a - tt 2 (a', r)d a ,)8(a - a'), 

where 5' (a — a') = -^5(a — a'). It should be noted that each of the above Poisson brackets 
leads to a set of terms at different points a and a' multiplied by -§^5{a — a') or -^j5{a — a') 
which equals to —J^5(a — a'). However, since these terms have only non vanishing value 
when a' approaches to a, one can consider all of them at the same point. Then they add 
up to give the above results. The algebra (fT3"j) shows that \I/ 1 (cr, r) and \J/ 2 (cr, r) are first 
class constraints. Moreover, from (jSJ) we see that: 

{V 1 ,H} = iVa'tt 1 + TV/* 2 + iA^tt 2 ' ^ , , 

{tt 2 , H} = iV/tt 1 + A^'tt 2 + iVxtt 1 + |iV 2 tt 2 ' « ^ 
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This shows that the consistency of \I/ 1 (cr, r) and \I /2 (cr, r) does not give any new constraint, 
and we are left with the five first class constraints given in ([T2"j) . 



In this way we have derived three constraint chains ( ^ j , ^ ^ 2 j and ( tt 3 ) 



in 



the terminology of reference [9]. In fact, the chain relation {it 01 ,!!} = holds for all of 
the chains. However the first two chains are correlated, since the Poisson bracket of the 
last element of each chain with the Hamiltonian contains the other constraint. This means 
that it is not possible to construct closed algebra within each chain. The last chain contains 
just one element and is not correlated to other chains, since it commutes with all of them 
as well as with Hamiltonian. 

As in ordinary Polyakov string one can show that 7r 3 generates the Weyl symmetry of 
the model which affects only the components of the world-sheet metric. In terms of the 
variables N a we have — > N 3 + e under Weyl transformation. On the other hand the con- 

/ tt 1 \ / 7T 2 \ 

straint chains I ^ 1 1,1 ^ 2 J can be shown that generate the effect of reparametrization 
invariance on the metric variables Ni and N 2 as well as the variables X a . 



3 Gauge fixing 

We began the theory with 14 field variables in the phase space, i.e. X a , N a and their 
corresponding momentum fields p a and ir a . Then we derived 5 first class constraints given 
in ffl2|) . As is well known from Dirac theory the first class constraints are generators of 
gauge transformations [10J. One needs to consider additional conditions to fix the gauges. 
These "gauge fixing conditions" are functions of phase space variables which should vanish 
to fix the gauges. The gauge fixing conditions should fulfill two conditions. First, they 
should constitute a system of second class constraints when added to the original first 
class constraints of the system. This condition is necessary to fix the values of variables 
which vary under the action of gauge generators [12]. Second, they should have a closed 
algebra under the consistency conditions, i.e. under the successive Poisson brackets with 
the Hamiltonian. 

For a " complete gauge fixing" the number of independent gauge fixing conditions should 
be equal to the number of first class constraints |13j . In this way, we should suggest 5 gauge 
fixing conditions to fix the gauges generated by the constraints given in ffl2l) . and reach 
a "reduced phase space" of 4 field variables. Since the momenta 7r Q are generators of 
transformations in N a , we fix the corresponding gauge by choosing the values of N a as 
N± ~ 1, ~ and N$ & 1. These values are chosen such that = rjij. In this way we 
have so far introduced three gauge fixing conditions 

ft x = Nx - 1, 

fi 2 = N 2 , (15) 

n 3 = N 3 -i. 

It can easily seen that the system of 6 constraints n a and Q a are second class. The consis- 
tency of Q a 's by the use of total Hamiltonian fllOp determines the lagrange multipliers X a 
to be zero and does not give any new constraint. This makes us sure that the two criterions 
of a good gauge mentioned above are satisfied. In fact, by the above gauge fixing three 
degrees of freedom N a are removed completely from the theory. This gauge has fixed the 
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Weyl symmetry as well as the effect of the reparametrization on the metric variables Aq and 
iV 2 . On the other hand, we are still left with the remaining gauges generated by ^f 1 and \I/ 2 
which generate the effect of reparametrization on the variables X a . In fact, since we have 
fixed the gauge from the middle of the constraint chains, the gauge is fixed partially in the 
language of reference [13]. In partial gauge fixing the Lagrange multipliers are determined 
while the variations generated by some of the gauge generators are not fixed. 

To fix the effect of the parametrization of the world-sheet on X a 's, as in so many 
models in string theory we need to determine some definite combinations of fields as the 
time variable in target space. Taking a look on the form of the constraints ^ and \I/ 2 
in OH]) shows that the choice u = pr is more economical in the sense that simplifies the 
constraints better. Here p is a parameter with dimension of (length) -1 . We recall that all 
of the dynamical variables in the action are dimensionless. Hence, we consider the gauge 
fixing condition 

Q4 = U — [IT. (16) 

To fulfill the second criterion of a good gauge we choose the last gauge fixing condition as 

= {^,H T } + ^ (17) 
~ Pv - 2/i 

This new constraint should also be valid during the time. Since 

No 
iVi 

the chosen gauges are consistent and make a closed algebra with the Hamiltonian. It is also 
clear that Q4 and Q5 make a second class system with ^Z 1 and \l/ 2 . Imposing strongly the 
constraints (TIB"]) and (fT7|) on the system, simplifies the constraints \Eq and ^2 as 

vj/ 1 -> fr 1 = Ipf + ±eijfj,aiPj + eijUra'iPj + (1 + p 2 r 2 )af + pp u - bp 2 + \p 2 af + p 2 Ta i a' i , 
\£ 2 _> \J 2 = a '.pi + 2pv', 

(19) 

This shows that p u and v can be derived on the constraint surface, i.e. from identities 
^1 = and ^2 = 0, in terms of the physical variables di and p^. In this way the reduced 
phase space is just the four dimensional space of (ai,pi) whose original Poisson brackets 
serve as the Dirac brackets in the remaining physical space. The terms pp u and p?b in the 
expressions of ^1 have nothing to do with the dynamics of («i,Pi) and can be dropped. The 
parameter b has in fact no important role in the theory and only shifts the spectrum of the 
energy with a constant value. 

As in reference [8] we consider the dimensionless quantity pi as a small parameter which 
should be considered only in the first order. Therefore, in all of the foregoing calculations 
we keep only linear terms with respect to p, assuming that I is finite. Therefore, the 
Hamiltonian (jSJ) in the reduced phase space can be written in terms of the Hamiltonian 
density: 

ri G F = -P 2 + ^pd^ + eijpra'tPj + af. (20) 

Since B(X) in ([2]) is linear with respect to u one may think of p as the order of magnitude 
of the 5-field. This assumption is equivalent to considering the effect of the 5-field only 
up to the first order. 
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4 Boundary conditions as constraints 



From now on we forget about the original theory and suppose we are given a theory with 
two degrees of freedom a, and the corresponding momenta whose dynamics is given by the 
final Hamiltonian fTSOj) . We make a change of variables from (cii,pi) to (Ai = e^aj, Pi = Pi). 
Then the the fundamental Poisson brackets which is the same as the final Dirac bracket of 
the original theory read 

{A(^,r),Pj(a',T)} = e lj 5(a - a'), 

{A(a,r),^K,r)} = {P^r), P>' »} = ^ 
The Hamiltonian equation of motion for the remaining fields, can be written as 

Ai = 2 € ij(Pj ~ 2/ItA'j — flAj) . . 

P i = -e ij (frP j - f JiTP , j + 2A'j) 1 } 

The only things that should be brought from the original theory are the boundary 
conditions. Using the original action (J4]) the boundary condition after gauge fixing emerge 
in terms of phase space variables as: 

= nrPi - 2A' i = at a = 0,1 (23) 

We have shown in the appendix that the boundary condition (123]) can also be derived from 
the parallel approach as the equations of motion of the end points in the discretized version. 

As mentioned in the introduction we do not want to find the general solution of the 
dynamical equations of motion. On the other hand, we are interested to follow the dynamics 
of the boundary conditions which means investigating the consistency of primary constraints 
$ J - 1 ' ) (cr)| (J= o and $^ (c)!^. Using the gauge fixed Hamiltonian of the previous section (|20|) 
the total Hamiltonian at this stage is 

Ht = J da[\PiPi ~ \^ p i - ^APi + AA] + Ai$fV)Uo + A^ (1 V)U*- (24) 
The consistency of primary constraints for instance at a = gives 

0= ^P ? ,-e 4J P; + ^4] CT=0 + A{ {$S\=o (25) 



Similar equations should be written at the end-point a = I. As discussed in details in [14] 
the first term in the LHS of Eq. (125j) has not the same order as the coefficient of (and A\) 
in the second term when regularizing the Dirac delta function. Therefore this condition can 
be fulfilled identically only if A^ 2 as well as the first term vanish simultaneously. In this way 
we have used the consistency conditions of the constraints for simultaneously determining 
the undetermined Lagrange multiplier and finding the next level of constraints as <&\ 2 \o) 

(2) 

and $j (/) where 

$f \a) = P t - e^Pj + iie^A' r (26) 

Then we should consider the consistency of second level constraints by using the Hamilto- 
nian ^ 

H = J dai-PiPi - -fiAiPi - ixrA\P t + A\A[] (27) 
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which is the same as the total Hamiltonian (J21]) after imposing A\ 2 = 0. This gives the 
third level of constraints. Subsequent levels of constraints can be derived in the same way. 
Using the relations: 

{A^,H} = i^(jf> " Mj n) - 2/^4^) + O(^) 

where Af ] = <9™ A* and P/ n) = d^Pi one can inductively show that the full set of constraints 
are $f°(0) « and « where 

.(2n+l) _ _, p (2n-l) , ,_p(2n) ~ ^ , (2n) Q , (2n+l) o\ 



f n+1) + (2n + 1 - \> 



$f"+ 2 ) = („ + i)^/ 2 ") - e .. Pf n+1) + (2n + l)/^ 2 "^ + 0(// 2 ) n = 0, 1, 2, 



(29) 

For practical calculations we write the constraints as ordinary functions in the bulk of the 
string and then integrate them with the use of 6(a) and 5 (a — I) respectively. 

Now we want to investigate whether the constraints are first or second class. For this 
reason one should calculate the Poisson brackets of the constraints. Since the constraints 
contain different orders of derivatives of Ai(a, r) and Pi(a, r), the Poisson brackets Cf- = 
{$f , } contain derivatives of orders k + k', k + k' — 1, etc, of the Dirac delta function, 
which are highly divergent and independent of each other. One way of treating the matrix 
of Poisson brackets is regularizing the delta functions as gaussian functions of width e and 
let £ — >■ after all. A tedious calculation gives 

C 2rn + l,2n+l = -^ £ _ 2(m+n+1) ^ + njH 2m+2n (0) - 2rH 2m+2n+1 (0)) + O(^) 
C 2rn + 2,2n + l = ^^(m+n+l)-!^..^^^ + Sij H 2m+2n+2 (0)) + 0(/l 2 ), (30) 
-i2n 



c 2rn + 2,2n+2 = ^^(^n+l)-!^^ (Q) + Q ^ 2) 



where H n (x) are Hermite polynomials. Similar expressions should be considered with H n {l) 
at the end-point a — I. The non vanishing elements on each row are located such that no 
vanishing linear combination of rows may be found. This means that the infinite dimen- 
sional matrix Cfj k is not singular and can in principle be inverted. Therefore, all of the 
constraints are second class. However, it is not practically possible to find the inverse of 
Cfj k . The problem is how we can find the Dirac brackets of the fields which need to have 



5 Reduced phase space 

As stated before, we seek for appropriate coordinates in which imposing the constraints 
( T2D"j) lead to omitting a set of canonical pairs. Here we have a difficult problem in which the 
ordinary Fourier expansion does not do this job. However, in the limit \x — > the boundary 
condition (123]) is the ordinary Neumann one and the Hamiltonian ( |2"7|) has a simple quadratic 
form in terms of coordinates and momenta. Hence, we need to write extended Fourier 
transformations for the fields Ai and Pj that include at most linear corrections with respect 
to the parameter /i and go to the ordinary Fourier transformation in the limit — > 0. Since 
fir and fia are the only dimensionless quantities that can be used for this correction, what 
can we do is correcting the Fourier coefficients by correction terms linear in r or a. The 
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linear term in r, however, is not needed at this stage, since it can be considered as part of 
the solution of the equations of motion. Adding all these points up together we suggest the 
following extended Fourier transformations for the fields 



1 



oo 



Ma, t) 



V2^ 



dk [(Ai(k, t) + fiaai(k, r)) cos ka + (Bi(k, r) + ua(3i(k, r)) sin ka] , 



— oo 



(31) 




V 7 ^ 




dk [(Cj(k, t) + fiajj(k, r)) cos ka + (Dj(k, t) + uaSj(k, r)) sin ka] . 



(32) 



In ordinary Fourier expansions the coefficients Ai(k, r), Bi(k, r), Ci(k, r) and Di(k, r) con- 
tain the same amount of data as the original fields Ai(a, r) and Pi(a, r). Comparing the 
expansions ( l3~Tj) and (I3"2"|) with ordinary Fourier expansions shows that we have used a du- 
plicated basis including sin's, cos's, a times sin's and a times cos's for expanding our fields. 
This basis is complete but its elements are not independent. Mathematically it is allowed 
to use a basis which is "larger than necessary". However, the essential point is that one 
should assume appropriate Poisson brackets among the extended Fourier modes such that 
the desired fundamental Poisson brackets ( 12T1) remain valid. In other words, we should 
tune their brackets in such a way that our physical phase space variables, which are half of 
the extended phase space variables, do obey the right Poisson brackets. Direct calculation 
shows that the following Poisson brackets lead to the standard Poisson algebra ( 121]) for the 
physical fields, 



All other Poisson brackets are assumed to vanish. Specially the modes j3i and 5i have 
vanishing Poisson brackets with all other variables in the extended Fourier space and so 
decouple from the theory. This means that we can put them away and write down the 
expansions only with linear terms in the cosine modes. We will see on the other hand 
that omitting the modes and 5i does not disturb our analysis of imposing the boundary 
conditions. We have, up to this point, 6 sets of real variables in the extended Fourier space 
which depend on real, continues and positive variable k. 

Now we want to impose the full set of constraints (|2"§|) on the fields. Using the expansions 
( l3"Tj) and ( 13"2"|) the constraints at the end-point a = lead to 



where Bi = uBi and Di = fiDi. Since these conditions should be satisfied for arbitrary 
values of n we have 



The difficulty arises here since the integer n, which shows the level of constraints, has 
appeared in the form of relations among the Fourier modes. This means that it is not 
possible to satisfy the constraints of all levels just by considering simple linear relations 



{^(fc,r),C i (fc',T)} = {Bi(k,r),Dj(k',r)} — 8ij8(k - k'), 
{a i (k,T),D j (k',T)} = {>y i (k,T) ) B j (V,T)} = 6 ij d k ,6{k-V). 



(33) 



dkk 2n fireijCj + 2ne lj A j + (An + 2)c^ + 2kB { + 0(fi 2 ) = 
JZz dkk 2n ~ l [(n + l^C, + (2to + l) 7i + kDi] + 0(fi 2 ) = 



(34) 



fireijCj + 2ne ij A j + (An + 2)a { + 2kBi = 0, 
(n + VjEijCj + (2to + l)/i7i + kDi = 0. 



(35) 



S 



among the Fourier modes of a given k as can be done in ordinary Dirichlet, Neumann, or 
even mixed boundary conditions [Bj. In fact, this phenomenon is the reason which makes 
the ordinary Fourier expansion inadequate for realizing the constraints. However, we have 
the opportunity of existence of extra variables in the extended phase space, which provides 
us additional tools for satisfying the constraints. In this way we are allowed to assume that 
the coefficients of n besides the terms independent of n in (133]) vanish. This gives 



a t = -\eijAj + C(/i 2 ) B, = ^-(A,- - tCj) + 0(fi , 



Hence the main fields Ai(a,r) and Pi(a,r) can be written in terms of two remaining sets 
of Fourier modes A^k, r) and Ci(k, r) as 



Ai(a,r) = -= / dk 

V Z7T 



(8ij — -fia6ij)Aj cos ka + 7£T e ij (A? — T Cj) sin ka 



Pi{a,r) = -= I dk 
v 2n 



{fy + 7 ) k L vSij)Cj cos ka + ^-Tpi sin ka 



(37) 



(3? 



As expected, the zeroth order (with respect to \i) of the Eqs. (]37[) and 038p is the expansion 
of a simple bosonic string with Neumann boundary condition at the end point a = 0. The 
linear term with respect to a in cosine modes as well as the sin term itself are appeared as 
the first order corrections. 

Next we should impose the constraints (1291) at the end-point a = I on the fields derived 
recently in Eqs. (J37l) and (J38l) . Hence we find 

/f dhk^i-lfinixeijCj + 2P(A i - i/i^)] sin(kl) + C(/x 2 ) = 0, 

(39) 

dkk^i-imSij - \^ae ij )C j - (2n + l)^A 3 ] sm(kl) + 0(^ 2 ) = 0. 

The above constraints are satisfied identically for kl = mn. However, for ^ / ™ there is 
no way for satisfying the constraints for arbitrary n except assuming that 

777 7T 

Mk,T) = Ci(k,T) = Q for (40) 

This leads to descritizing the Fourier modes. 

Before writing the final form of the fields in terms of the set of enumerable Fourier 
modes, care is needed to write the zero modes. The contributions due to cosine modes 
come out automatically by letting k — 0. However, contributions to zero modes originating 
from sine terms should be derived by taking the following limits: 

1 ~ 1 

lim Bi sin ka = -aeij(Aj(0, r) — rCj(0, r)), lim Di sin ka = — a€ijCj(0, r), (41) 

which follow from Eqs.( l36p . Adding these two contributions the zero mode part of the fields 
are so far as follows 

A°(a,T) = A°(T)-^are li q(T) 

P?(a,T) = -(e l] +va5 lj )q(r) ^ 

At this point we want to notice the reader to a global symmetry of the gauged fixed 
Lagrangian. If we turn off the B-field we would have an ordinary bosonic string in which 
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only the derivatives of the A-fields are present in the Lagrangian. This allows one to shift 
the fields by a constant amount without any change in the Lagrangian. When the B-field 
is on, Eq. ( l2"Uj) shows that the A-field itself is present in the gauged fixed Hamiltonian. 
However, the relevant term, i.e. the second term in Eq. (|20|) . is proportional to fx. This 
shows that the theory is symmetric, up to second order terms with respect to fx, under the 
following transformation 

M<r,T)^M<r,T) + fif(r) (43) 

where /(t) is an arbitrary function of time. This symmetry leads to an ambiguity in the 
zero mode of the A-field. Hence we should correct the first row of Eq. f!42p in the most 
general CcLSG clS follows 



4V, r) = A°(r) - -fxare^r) + fxl^A^r) + 6yCj(r)] 



(44) 



Note that fxl is the only relevant dimensionless quantity which is first order in fx. The 
unknown coefficients and bij should be determined upon suitable assumptions about 
the algebra of the fields. The best assumption seems to be keeping the standard algebra 
( 12 lj) in the bulk of the string and letting all changes in the algebra of the fields lay on the 
boundaries. If we make this choice the final form of the physical fields in terms of the set 
of discrete Fourier modes Afir) = AA^,t) and CT(r) = CA^t) are as follows 



Ada, t) 



j_ 



+ 



4V) - 

2 v^oo 
I Z^m=l 



fxr(a - ^q(r) - l^A^r) 



(A™(t) - ^ae rj Af(r)) cos ^ + ^ €ij (Af(r) - rCf(r)) sin : 



T 



^t))cos^P 



^C™(r) sin^ 

2n.7r IV/ I 



The normalization factor 



r is replaced by ^ j for oscillatory modes and 



(45) 



(46) 



for zero 



mode upon going from the continues parameter k to the discrete parameter m. B With this 
normalization the brackets of the discrete modes should also be given in terms of Kronecker 
delta as 

{A?,Cf } = 8 ij S mm ,, (47) 
{A?,Af} = {Cr,Cf} = 0. (48) 

In fact, the remaining canonical pairs A™ and Cf 1 as a small part of the original phase 
space are natural coordinates of the reduced phase space. On the other hand, a great part 
of the initial phase space variables are omitted due to the constraints. 

Remember that if one is able to omit the redundant variables due to all kinds of con- 
straints and write down the relevant fields in terms of final canonical coordinates of the 
reduced phase space, then there is no need to find the Dirac brackets. In other words, we 
pay the expense of using the Dirac brackets whenever it is not possible to find a canonical 



3 Since another length scale, i.e. exists in the model, one may suppose that the normalization 

factors should differ from the ordinary Fourier scries. However, it can be shown that such corrections only 
changes the observables by amounts of 0(fi 2 ) which is not important 
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basis to describe the reduced phase space. Hence, we will find the Dirac brackets of the 
original fields A4 (a, r) and Pi(cr, r) if we calculate their brackets by using the brackets (j4"7|) 
and flUD- 

Eq. (|46|) shows that the momentum-fields Pj(<x, r) just include the variables CJ™ and 
have vanishing brackets: 

{P l (a,r),P J (a',r)} = 0. (49) 
Straightforward calculations gives the brackets of coordinate and momentum fields as 

{Ai(a, T),Pj(a f , t)} = eij8 N (a, a'), (50) 

where 

S N {a, a') = 5{a - a') + 5(a + a'). 

Since both a and a' lie in the interval [0,/] their sum never vanishes. So the second 
delta function does not have any role and Eq. ( 150]) reduces to the usual form of Eq. 
(12 ip . However, since in the expansion of A-fields both variables A™ and C™ 1 are present, 
the interesting phenomenon appears in the bracket of coordinate fields at different points. 
Direct calculation gives 



{Ma, r), A 3 (a', r)} = -/xre^ I - 1 



2 ^ 1 . n-K , \ 

- > -sin— [cr + a) 
7T ^ n I J 

n=l / 



This result is similar to what derived in [6] for a string coupled to constant background 
B-field. The right hand side of Eq. ( |5T1) vanishes in the bulk of the string, i.e. when 
a or a' does not lie on the end points. It gives (-2) when a = a' = and (+2) when 
a — a' — I. However, as the B-field itself, the amount of non commutativity grows linearly 
with time. Our result here defers from reference [7j with a term proportional to /if 2 which 
is the same on both boundaries as well as in the bulk of the string. If, however, we add a 
term —}^\xt 2 e^C^-ir) to the zero mode part of the field A;(a, r) in Eq. ( I4"5j) . our result will 
coincide with reference [7]. This correction is allowed according to the global symmetry of 
Eq. (T4"5|) . This means that we have forgiven our previous assumption that the components 
of the A-field commute in the bulk of the string. With this assumption the resulted brackets 
can be summarized as follows 

{A i {<r,T),P j (o',T)}-- 

{Pi(*,T),P j (o',T)} = 

a ^0,1 or a' ^0,1 (52) 
{Afor), Mo'it)} = { (xre^l + £) a = a' = 

a = a' = I 

This shows that the fundamental characters of the A-fields and P-fields as coordinate and 
momentum fields are remained almost as before and the time dependent B-field leads to a 
time dependent non commutativity in the coordinate fields all over the string. 




6 Concluding remarks 

In this paper we gave a complete Hamiltonian treatment of the Nappi-Witten model (WZW 
model based on non semi simple gauge group) as an interesting and non trivial system in 
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which complicated boundary conditions make the physical subset of variables far from 
reaching. The initial dynamical variables in this model are 4 components of a bosonic 
string, X a = (ai, a 2 ,u,v), and the components of world-sheet metric. We used appropriate 
variables to find 3 primary and 2 secondary first class constraints. It can be shown that 
these constraints are generators of reparametrizations as well as Weyl transformations. 
Then we fixed the gauge such that the world-sheet metric is flat and u = ^lt where the 
small parameter \i is proportional to the strength of the B-field. In this way the components 
of the world-sheet metric and the variables u and v disappeared as the result of constraints 
and gauge fixing conditions. Hence, we derived a smaller theory with two coordinate fields 
ai and ci2 and their corresponding momentum fields. 

The most important part of the problem seems to be the boundary conditions which 
should be brought from the original theory. Considering the boundary condition as Dirac 
constraints and following their consistency, we found two infinite chains of second class 
constraints at the end-points which restricted the space of physical variables to a much 
smaller set. Due to complicated form of the boundary conditions, it is not an easy task 
to impose them on the space of the physical variables. In fact, with an ordinary Fourier 
expansion the constraints do not lead simply to omitting some Fourier modes as in Dirichlet 
or Neumann boundary conditions. 

To overcome this difficulty we extended the phase space to a larger one which is given by 
an extended Fourier expansion in which the Fourier modes are replaced by linear functions 
of the variables. In this basis the infinite set of constraints can be imposed more easily 
by using the arbitrariness due to extra variables. This results to disappearing of so many 
canonical pairs among the used extended Fourier basis and finally a set of discrete modes 
remain which act as the canonical coordinates of the reduced phase space. Then all physical 
objects including the original coordinate and momentum fields can be expanded in terms 
of these modes. 

Using these expansions we found that the commutation relations of the coordinate and 
momentum fields are almost as usual, except that the coordinate fields do not commute 
at the boundaries, with an amount proportional to time and/or B-field but with opposite 
signs at two boundaries. We showed that it is allowed to insert a term which gives non 
commutativity proportional to r 2 throughout the string. This correction may make our 
results consistent with those of reference [7] in which the authors have given iterative 
solutions for the equations of motion. 

We think that our method here has two main advantages in two different areas. First, 
we do not solve the equation of motion. Therefore, in our final result the time dependence 
of remaining modes are not specified. However, this time dependence is not essential for 
quantization of the model. If needed, one can use the Hamiltonian written in terms of the 
final modes and then derive their time dependence. In fact, our main objective is that for 
quantizing a theory, i.e. investigating the algebraic structure of the observables, it is not 
needed to follow the full dynamics of the system; it is just enough to study the dynamics of 
constraints. As a matter of fact, for simple models it may seem more simple and economic 
to solve the equations of motion and then quantize the theory, since this procedure contains 
the dynamics of the constraints within itself. But this may not be the case for a complicated 
model such as the model considered in this paper. 

The next advantage is in the context of constraint systems. As we see in the literature 
[U [H] the main difficulty in considering the infinite set of constraints due to boundary 
conditions is deriving the Dirac brackets. In this paper, as in our previous work [6] we 
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showed that if one is able to find a set of canonical variables describing the reduced phase 
space, then there is naturally no need to calculate the Dirac brackets. In fact, this was the 
main brilliant idea of Dirac [TT], who gave his famous formula of Dirac brackets in such a 
way that it is equivalent to calculating the Poisson brackets only in the space of canonical 
variables describing the reduced phase space. 
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